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ABSTRACT 

Context. Planning of the Square Kilometre Array (SKA) requires simulations of the expected sky emission at arcsec angular res- 
olution to evaluate its scientific potential, to constrain its technical realization in the best possible way and to guide the observing 
strategy. 

Aims. We simulate high-resolution total intensity, polarization and rotation measure (RM) maps of selected fields based on our recent 
global 3D-model of Galactic emission. 

Methods. Simulations of diffuse Galactic emission are conducted using the hammurabi code modified for arcsec angular resolution 
patches towards various Galactic directions. The random magnetic field components are set to follow a Kolmogorov-like power law 
spectrum. 

Results. We present maps for various Galactic longitudes and latitudes at 1.4 GHz, which is the frequency where deep SKA sur- 
veys are proposed. The maps are about 1!5 in size and have an angular resolution of about 1"6. We analyse the maps in terms of 
their probability density functions (PDFs) and structure functions. Total intensity emission is more smooth in the plane than at high 
latitudes due to the different contributions from the regular and random magnetic field. The high latitude fields show more extended 
polarized emission and RM structures than those in the plane, where patchy emission structures on very small scales dominate. The 
RM PDFs in the plane are close to Gaussians, but clearly deviate from that at high latitudes. The RM structure functions show smaller 
amplitudes and steeper slopes towards high latitudes. These results emerge from the fact that much more turbulent cells are passed 
through by the line-of-sights in the plane. Although the simulated random magnetic field components distribute in 3D, the magnetic 
field spectrum extracted from the structure functions of RMs conforms to 2D in the plane and approaches 3D at high latitudes. This 
is partly related to the outer scale of the turbulent magnetic field, but mainly to the different lengths of the line-of-sights. 
Conclusions. The significant scatter of the simulated RM distributions agrees with the large scatter of observed RMs of pulsars and 
extragalactic sources in the Galactic plane and also at high Galactic latitudes. A very dense grid of RMs from extragalactic sources 
is required to trace and separate Galactic RM foreground fluctuations. Even at high latitudes total intensity and polarized emission is 
highly structured, which will contaminate sensitive high-resolution extragalactic observations with the SKA and have to be separated 
in an appropriate way. 

Key words, polarization - Radiation mechanisms: non-thermal - ISM: magnetic field - ISM: structure 

1 . Introduction used to find their distances based on a thermal gas distribution 

model. Recent pulsar results indi cate a num ber of ma gnetic field 

Intensive radio emission from the Galaxy originates from ion- reyersals wkhin ^ Gal ( |Han et alf 2006: NoutsosIFID 

ized thermal gas, cosmic rays and magnetic fields. All-sky sur- j^g^ which h difficult tQ understand in terms of a si le axi . 

H ; veys yield an overview of the Galactic structures. However, to etric or bi _ symme tric magnetic field configuration. A more 

reveal distributions of the Galactic emission components, mod- denfje RM id h needed tQ deriye a ^ l cale ttem> 

elling is required because of our position inside the Galaxy. becaufje lications arise from local magnet ic field distur- 

Galactic emission is a severe foreground contamination for ob- bancef , Qr excessive thermal gas contributions t0 the observed 

serrations of fluctuations of the cosmic microwave background , sar DMs and RM§ (Mkra gt aJ 2003) Such effects must fee 

radiation and of all kinds of distant, faint extragalactic objects, se ated to derive a reliable Yalge . scale magnetic field config- 

which will be investigated by the Square Kilometre Array (SKA) uration Moreoveri t0 deC ompose the observed RM information 

with unprecedented sensitivity. Therefore it is crucial to prop- mtQ a netic field contribution , a det ailed knowledge of the 

erly account for Galactic emission in all directions of the sky. themal distribution is required . The currently widely ac- 

Unfortunately, most of the properties of the aforementioned con- d distribution of Galactic diffuse thermal gas, the NE2001 

stituents of the Galaxy from large scales to small scales are in- model> afj weU af , ^ fim factQr m odelled fr om pulsar and 

sufficiently constrained up to now. Hff observations was prese nted by Cordes & Lazio ( 2002) and 

The large-scale Galactic magnetic field is of particular in- B erkhuijsenetalJ dS), respectively. Galactic magnetic field 

terest. It has been intensively studied by using rotation mea- models based on pulsar m observations have to agree with the 

sures (RMs) of pulsars (e.g.| Han et al.|2006HNoutsos et al.|2008|) observed rm s of extragalactic sources shining through the entire 

and extragalactic sources (EGSs) in the Galactic plane (e.g. Galaxy To establish a firm picture of the i arge . SC ale magnetic 

| Brownetal. || 2003i | 2007|) . Pulsar dispersion measures (DMs) are fieM needs tQ acco(mt for ^ me effects aboye 

Send offprint requests to: Wolfgang Reich Recently, ISun et alj (I2008I) have derived a Galactic 3D- 

e-mail: wreich@mpifr-bonn.mpg.de emission model which reproduces all-sky surveys over a wide 
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frequency range and also the all-sky RM distribution of EGSs. 
The simulated maps have an angular resolution of about 15', 
which is suitable to compare them with all-sky maps having 
lower angular resolution. Based on these global results for the 
Galaxy, simulated arcsec angular resolution maps are the next 
step in order to get an idea what might be observed with the SKA 
and other future high-resolution radio telescopes. It is important 
to understand the origin of the observed structures at high angu- 
lar resolution to estimate the contamination of extragalactic ob- 
servations by Galactic small-scale structures as the foreground. 

One of the key science projects for the SKA is "Cosmic 
Magnetism", where RMs of about 2 x 10 7 polarized EGSs 
are proposed to be o bserved in a region of 10 4 deg 2 
(Beck&Gaensler 2004). The survey is proposed for 1.4 GHz 
with a resolution of about 1". Simulations of the expected num- 
ber density of faint sources and their p olarization propertie s 
from SKA observations were made by e.g. lWilman et al.l (12008). 
A feasibility study on the decomposition of the source RMs 
from the Galactic foreground RMs and its removal need high- 
resolution Galactic RM maps. This is the basic motivation for 
extending our simulations to selected Galactic patches with arc- 
sec angular resolution. 

Fluctuations of the components of the interstellar medium 
(ISM) dominate the emission structures at high angular resolu- 
tion. The thermal electron density distribution has been found 
to follow a Ko lmogorov-like power l aw from about 10 8 m to 
about 10 18 m (lArmstrong et al.lll995l) . The fluctuations of the 
Galactic magnetic field are mainly studied based on the anal- 
ysis of the structure function of RMs, where an additional ef- 
fort is needed to de couple the magnetic field and the electron 
density fluctuations (Mi nter & Spangler]|1996l) . The outer scale 
of turbulence was reported to be in the ran ge of about 100 pc 
( lArmstrong et al.ll 19951: Haverkorn et al . 2008) down to a few pc 
dMinter & Spanglerll 19961: IHaverkorn et al.ll2008l) . 

We summarize our global Galactic 3D-modelling in Sect. 2. 
The high angular resolution simulation method and the treatment 
of small-scale magnetic fields are described in Sect. 3. Results 
and analysis of the simulated RM maps are presented in Sect. 4, 
including a discussion of the consequences of the simulated RM 
fluctuations. In Sect. 5 we discuss the results of the simulated 
total intensity and polarization maps. Some general conclusions 
are given in Sect. 6. 



2. A new Galactic 3D-emission model 



I Sun et all d2008l) presented a new Galactic 3D-emission model, 
which is in agreement with a wide range of observations: 

- Galactic synchrotron radi ation being best repr esented by the 
408 MHz all-sky survey dHaslam et al.ll 19821) in total inten- 
sity and b y the WMAP 22 .8 GHz all-sky survey in polarized 
intensity (Pag e et al.ll2007l) . 

- Optically thin thermal emis sion by the WMAP 22.8 GHz 
free-free emission template (Hins haw et alj |2007). and opti- 
cally thick thermal emission causing absorption, which re- 
flects in a spectral flattening at low frequencies along the 
Galactic plane. 

- The all-sk y 1420 MHz polariza tion survey combined from 
northern ( Wolleben et al.l 120061) and southern sky data 
dTestori et al.l2008h constraining depolarization properties of 
the interstellar medium. 

- RM data from EGSs along the Galactic plane taken from the 
Canadia n (CGPS) and the South ern Galactic Plane Survey 
(SGPS) dBrown et alj2003ll2007l) . and EGSs out of the plane 



from lHan et al.l d!999l) . including preliminary RMs from the 
new Effelsberg L-band survey (Han et al., in prep.), and 
RMs from the CGPS high latitude extension (Brown et al., 
in prep.). 

The Galactic 3D-model is based on the 3D-distribution of 
thermal electro ns in the Galaxy (NE20 01 model) derived from 
pulsar DMs by ICordes & Laziol d2002l) . To successfully model 
the observed thermal emission for the optically thin case and 
absorption effects from the optically thick gas, the filling fac- 
tor of thermal gas needs to be taken into ac c ount, whose spatial 
variation was derived by Berkhuijsen et al. (2006). The regular 
Galactic disk magn etic field model wa s cons trained by extra- 
galactic RM data of iBrown et all d2003l 120071) . Axi-symmetric 
(ASS) and bi-symmetric (BSS) fields with one field reversal in- 
side the solar circle agree with extragalactic RM data in the 
plane, while the BSS field configuration fails to represent the 
RM gradient with Galactic latitude as measured from the high- 
latitude extensio n of the CGPS by B rown et al. (in prep.). The 
408 MHz survey dHaslam et al J 1 982l) and th e WMAP 22.8 GHz 
linear polarization survey dPage et alj|2007l) were used to con- 
strain the random magnetic field component, the cosmic-ray 
electron distribution, and the known local synchrotron excess. 
Finally, the depolarization effects of the 1420 MH z polariza- 
tion survey dWolleben et al.l2006tlTestori et al.l2008l) were mod- 
elled, which requires to modify the thermal filling factor by 
adding an extra coupling term of the electron density and the 
random magnetic field. 

The quantitative model parameters of the three constituents 
of the Galaxy are those of the NE2001 thermal electron density 
model, an ASS magnetic field configuration with a local strength 
of 2 fiG, asymmetric toroidal halo fields with respect to the plane 
with a maximum strength of 10 fiG, isotropic and homogeneous 
random magnetic fields with a strength of 3 fiG, and a cosmic- 
ray electron distribution trunc ated at 1 kpc abov e and below the 
plane. The details are given bv lSun et all d2008l) . 

Although this Galactic 3D-model is in agreement with the 
observed all-sky maps, some of its parameters seem unrealis- 
tic. To fit the observed high-latitude RMs, a strong regular halo 
magnetic field is required, because the scale height of the ther- 
mal electrons of the NE2001 model is just about 1 kpc. The 
strong halo field in turn requires a cosmic-ray electron distri- 
bution limit of about 1 kpc to avoid excessive high-latitude syn- 
chrotron emission. These parameters are un realistic for both the 
halo field strength dMoss & Sokoloflll2008l) and the distribution 
of cosmic-ray electrons. ISun et al.l d2008l) note that by increas- 
ing the scale height of the thermal electron density to about 
2 kpc the halo magnetic field strength reduces to about 2 fiG 
and the tr uncation of the cosm ic -ray distribution turns obsolete. 
Recently, Gae nsler et alj d2008l) critically reanalyzed the thermal 
electron scale height and quot e a revised value of 1.8 kpc, very 
close to the proposed 2 kpc bv lSun etal (2008). This makes the 
Galactic properties similar to those observed for nearby galax- 
ies. It should be noted that the simulated maps are not affected 
by these modifications of parameters. 



3. High-resolution simulations 

3.1. The general simulation concept 

The simulatio ns presented in this paper make use of the Ham- 
murabi code (Waelkens 2005;: IWaelkens et alj|2009j). where the 
sky is pixelized according to the HEALPix scheme dGorski et ail 
2005). Each pixel contains information from a cone shaped vol- 
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ume centered at the observer. Each surface element is quadri- 
lateral and covers the same area. The number of pixels N p ix of 
a sphere are calculated by N plx = 12A^ IDE , where A^ide is an 
integer power of 2. Therefore d0, the linear size or the angular 
resolution, can be calculated as 

d®=J^. (1) 

V 7T AfclDE 

For each pixel of the map, the Stokes parameters (/, U and 
Q) and the RM are integrated in the radia l direction fro m the 
observer up to a maximal distance r max dSun et alj |2008). The 
modelled emission distributions in the Galaxy are represented 
in cylindrical coordinates with the origin at the Galactic cen- 
ter. The distance to the Sun is taken as 8.5 kpc. Details of the 
regular magnetic field component, the cosmic -ray density and 
the thermal gas distributions in th e disk and halo of the Galaxy 
are described bv lSun et alj (|2008). The integration process sam- 
ples these distributions and inherently accounts for depolariza- 
tion. To accurately calculate the integral, each pixel is separated 
into many volume units with equal radial interval Ar. Thus the 
size of a unit can be approximated as rdOxrdOxAr with r being 
the distance of the volume unit to the observer. This process is 
repeated for all pixels, which are finally combined to an all-sky 
map. 

3.2. Zoom-in for specified patches 

According to Eq. (Q]) the possible resolutions depend on A^ide- 
To achieve arcsec resolution and to optimize the amount of com- 
putation time, we use A^ide of 131072, which corresponds to a 
resolution of about l'/6. The pixel number for an all-sky map 
then amounts to about 2 x 10 , which is far beyond any com- 
puter memory capacity today. Therefore high-resolution simula- 
tions are just possible for a small patch of sky. Thus we adapted 

the hammurabi code to be able to simulate parts of the sky. 

In the hammurabi code the NEST ordering dGorski et alJ 

120051) is used, which facilitates the zoom-in of specified patches 
as it is explained below. For A^side = 1 the entire sky is split into 
12 patches as shown in Fig.Q] They are numbered from to 11. 
For A^side = 2, each pixel is further split into four child pixels, 
also shown in Fig. Q] for patch "5". This patch is centered at lon- 
gitude of 90° in the plane with a diagonal size of about 90° x90°. 
Subsequently for A^ide increasing by 2 each pixel is further split 
into four pixels. In this scheme one can determine a target patch 
and make proper splits depending on the resolution. To reach a 
smaller patch size and higher angular resolution a much larger 
value for A^ide is required. 




Fig.l. The HEALPix pixel scheme dGorski et al l 120051) for 
Galactic coordinates used in this paper. 



3.3. Realization of the random magnetic field component 

It is generally accepted that the power spectrum of the ran- 
dom magnetic field component should follow a power law as 
P(k) oc k a . Here k - 2n/l is the wave vector and / is the spatial 
scale. For a Kolmogorov-like turbulence spectrum the spectral 
index is a = -11/3, which is used in all present simulations. 
Note that independent fluctuations of the electron density are not 
included in the simulations. Howver, fluctuations of the electron 
density are also expecte d to follow a Kolmogorov-like spectrum 
dArmstrong et al.iri995l) . The coupling of electron density and 
random magnetic fields we introduce d in our all-sky si mulations 
to model the observed depolarization dSun et al.l ?2008) results in 
Kolmogorov-like fluctuations for both components. 

Random magn etic field c omponents realized in the simula- 
tions presented bv lSun e t al. (2008) are generated in situ for ev- 
ery volume unit by Gaussian random number^]. To obtain the 
random magnetic field components following a given power-law 
spectrum, an inverse Fourier transformation must be performed. 
All available fast Fourier transformation (FFT) algorithms are 
based on a Cartesian coordinate system. We therefore set-up a 
cubic box with a size of l) and divide it into A^ 3 small cubes 
with a common volume of A 3 . Within the box we construct a 
Cartesian coordinate system and perform the FFT to derive the 
random field strength and direction in each small cube. Here L 
corresponds to the outer scale of the turbulence and A corre- 
sponds to the inner scale. 

An important constraint for A is that A < rdO. This ensures 
that several neighbouring line-of-sights do not pass the same 
small cubes. Otherwise spurious stripes are seen in the maps. 
For r of about 1 kpc and d© of 1'.'6 the size of a small cube 
should not be larger than about 0.01 pc. If our Galaxy is treated 
as one box, meaning L — 40 kpc, the total number of cubes is 
about 6.4 x 10 19 . Again there is no way in processing this huge 
amount of cubes by current computing facilities. 

In our simulation, we devise a scheme to generate the ran- 
dom magnetic field components in an adequate way. We first 
conduct a FFT to create a magnetic field distribution in a box. 
We then rotate this box to acquire new boxes, which are glued 
together to fill the entire Galaxy. The box rotation angles are in- 
teger multiples of n/2, which avoids cutting of boxes. Through 
this scheme the turbulence status of the Galaxy is realized. This 
process is detailed in AppendixlAl 

Actually the size of the box cannot be too small, otherwise 
there are too many boxes to handle. On the other hand L cannot 
be too large in order to get a reasonable number of cubes of about 
0.01 pc in size. As a compromise, we take a size of 10 pc for the 
box and of 0.0125 pc for the cubes. Although 0.0125 pc is much 
larger than the inner scale of the elec tron density fluctuatio ns, 
which are of the order of 10 8 cm (e.g. Armstron g et al.|[T995l) . it 
is sufficient for our study in regard to the current angular resolu- 
tion. In total there are about 5 x 10 8 cubes, requiring a memory 
of about 12 giga-byte for the random magnetic field components 
when stored as double precision numbers. It is thus obvious that 
the number of cubes is a compromise between high resolution, 
which requires small cube sizes, and hence a very large number 
of cubes and the available computer memory, which is thus able 
to handle a limited number of cubes. 

In Fig. |2] we show the 2D-distribution of random magnetic 
field components for a slice parallel to one surface of a 10 pc 
box. For comparison, we show the random magnetic fields gen- 

1 Throughout all the simulations we use the random number genera- 
tor from the GNU Scientific Library. The algorithm is MT19937, which 
has a Mersenne prime period of 2 19937 - 1. 
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Fig. 2. 2D random magnetic field distribution in a slice parallel to one surface of a 10 pc box. The upper left panel shows magnetic 
field components following a power law index of -11/3, while the upper right panel shows random field components realized by 
random Gaussian numbers corresponding to a spectral index of zero. Their corresponding power spectra are shown in the lower 
panels. 



erated in situ in each grid of the same pla ne as random Gau ssian 
number, which is similar to that used bv lSun et al.l ([2008). The 
power spectra of the magnetic fields are also shown in Fig. [2] It 
is apparent that the distribution of random magnetic fields fol- 
lowing a power-law spectrum shows coherent structures of var- 
ious extent, whereas the Gaussian random magnetic field distri- 
bution is featureless. The spectral index of ff2D = -2.83 ~ a + 1 
for Kolomogorov magnetic field distributions and <i2D = 
for Gaussian field distributions agree with the expectations and 
prove that the random magnetic field components are properly 
distributed in the box. 

4. Modelling results on rotation measures 

Except for the random magnetic fields the simulations presented 
below use the s ame Galactic 3D-emission models as derived by 
ISun et alJ d2008l) . where the "ASS+RING" model for the disk 
magnetic field is selected here. The random magnetic fields are 
introduced as described above and follow a Kolmogorov-like 
spectrum. The strength of the random fields is scaled to 3 fiG 



which is the same as that in ISun et all (2008). All simulations 
are made for 1 .4 GHz observations, as this is the proposed SKA 
RM survey frequency. As mentioned above, the simulated patch 
size is limited by available computer memory, so that just rela- 
tively small patches of about 1°5 in size are simulated at a reso- 
lution of 1'.'6. This is compatible with the expected SKA angular 
resolution. 

Simulations are presented for a handful of patches only, 
because each run takes very long computation time. We se- 
lect the first patch around (7, b) = (135°, +40°) which is in 
the d irection of the DRAO Planck Deep Fields ( Tavlo ret al.l 
2007), where soon high-resolution continuum and polariza- 
tion data will be available. This patch is named MN stand- 
ing for Medium-latitude-North. More patches with roughly the 
same longitude, but with different latitudes, are further se- 
lected. These are labeled HN (High-latitude-North), PNAC1 
(Plane-North- Anti-Centre 1), PS AC (Plane-South-Anti-Centre), 
MS (Medium-latitude-South) and HS (High-latitude-South). To 
compare with PNAC1, three additional patches in the plane, 
PNAC2 (Plane-North- Anti-Centre 2), PNC1 and PNC2 (Plane- 
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North-Centre 1 and 2) are included. The centre coordinates of 
these patches are listed in Table Q] Note that we cannot set the 
patch coordinates to exactly the same longitudes or latitudes be- 
cause of HEALPix constraints (Gorski et al. 2005). In Table Q] 
also the maximal integral lengths or the length of line-of-sights 
(/max) and the integral steps (Ar) are listed. The integral length 
is determined by taking into account the cosmic-ray scale height 
of about 1 kpc and the disk cutoff rad ius of 17.5 kpc for th e 
electron density in the NE2001 model dCordes & Lazioll2002T) . 
Integrating further out results in negligible contributions. For the 
patches with large integral length, the step size is larger than 
for those with small length to limit the computation time. The 
selected patches allow to study variations of emission and RM 
structures related to turbulent magnetic field components ver- 
sus Galactic longitudes and latitudes. Variations of the structure 
functions for RMs from E GSs as a function of Galactic direction 
have already been noted bv lSun &Hanl J2004) for angular scales 
of several degrees. 



Table 1. Simulated patches. 


name 


1(1 


b(°) 


r raa x (kpc) 


Ar (pc) 


HN 


138.33 


70.16 


3 


1.0 


MN 


136.48 


44.20 


3.5 


1.0 


PNAC1 


130.07 


1.19 


10 


2.5 


PSAC 


130.07 


-1.19 


10 


2.5 


MS 


136.48 


-44.20 


3.5 


1.0 


HS 


138.33 


-70.16 


3 


1.0 


PNC1 


59.77 


1.19 


30 


3.8 


PNC2 


300.23 


1.19 


30 


3.8 


PNAC2 


229.92 


1.19 


10 


2.5 



All simulations were made by using the computer cluster 
of the Max-Planck-Institut fur Radioastronomie operated by the 
VLB I group. 



4.2. Statistical description of the simulations 
The rotation measure, RM, is defined as 



RM = K 



Jo 



(flu + 6|,)dr, 



(2) 



where K is a constant, n e is the thermal electron density, By and 
b\\ are the regular and random magnetic field components pro- 
jected along the line-of-sight, respectively, and r max is the maxi- 
mal distance along the line-of-sight from the observer. One ma- 
jor question is to what extent the properties of B and b can be 
recovered from the RM maps. 



4.2.1. The probability distribution function 

We use the probability distribution function (PD F) to describe 
the str ucture properties of a map as discussed by iMizeva et al.l 
(I2007I) . For the RM distribution PDF (RM) x ARM quantifies the 
relative number of pixels falling between RM and RM + ARM. 
A Gaussian shape of the PDF indicates that no extended features 
are present in the map. To measure the deviation from a Gaussian 
distributi on, two more para meters, skewness and kurtosis, are 
used (e.g. lKowal et~aT1l2007l) . 



skewness = 



kurtosis 



1 y(rm,-rm 

1 y ( RM, - RM 



(3) 



4.1. Overview of the results 



where RM and <x RM are the average and standard variance of 
the RM maps. They are third and fourth order statistical mo- 
ments, respectively. Large absolute values of skewness or kur- 
tosis indicate that either the tail or the centre of the Gaussian 
is deformed. Physical ly these parameter s are closely related to 
turbulence properties (Kowal et al. 2007)). 

The RM PDFs for all simulated patched are shown in Figs. [4] 
[6] [7] and [8] (upper panels). Their average values, together with 
variance, skewness and kurtosis are listed in Table [2] 



As described by lSun etal the hammurabi simulations re- 

sult in Stokes /, U, and Q maps at the specified frequency as 
well as the RM maps. Subsequently, the polarization angle (PA) 
and polarized intensity (PI) maps are calculated according to 
PA = ±atan| and PI = ^U 2 + Q 2 . 

In this paper, we will not discuss the separation process of 
total and polarized intensities of EGSs from the foreground, 
which requires adding polarized EGSs and noise to the simu- 
lated maps. The simulated RMs represent the contribution of the 
whole Galaxy to the RM from an EGS in case the source size 
does not exceed the pixel size of the simulated map. Thus, the 
RM maps are directly related to observations. Below we focus 
on the analysis of RM maps, which are directly connected to the 
planned SKA RM survey. All simulated RM maps are shown in 
Fig. [3] The other maps are shown in Figs. IB.lllB.21IB.3IIB.4l and 
lB31 in AppendixE 

The RM maps clearly show structures ranging from small 
to large scales. The variation of the structures versus latitudes 
is also evident. Below some statistical methods are applied to 
analyse these structures. 



Table 2. Statistical parameters for the simulated RM patches. 
The patch designations are listed in Table 1. Note that for the 
patches HN, MN and MS more realizations were made with de- 
tails given in the text. 



name 


average 
(rad irT 2 ) 


variance 
(rad nT 2 ) 


skewness 


kurtosis 


HN 


2.50 


12.73 


0.0857 


-0.2555 


HN (Gaussian) 


2.14 


8.11 


-0.0007 


0.0008 


MN 


-1.90 


15.12 


0.0347 


-0.2287 


MN (Gaussian) 


1.41 


10.93 


-0.0007 


0.0009 


PNAC1 


-85.82 


63.40 


0.0654 


0.2124 


PSAC 


-103.13 


60.86 


-0.0482 


-0.0274 


MS (seedl) 


-27.58 


18.81 


-0.1464 


-0.6365 


MS (seed2) 


-19.78 


16.74 


0.1338 


-0.2968 


MS (seed3) 


-24.86 


16.35 


-0.1847 


0.3896 


HS 


1.66 


17.15 


-0.1241 


-0.6353 


PNC1 


-89.03 


93.65 


-0.0018 


0.0646 


PNC2 


286.53 


155.71 


-0.0045 


0.0528 


PNAC2 


159.80 


50.80 


0.0217 


-0.0429 






HS; RM (138,3!. 70 r l6) 
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HE5-; RH (3*b,M. l r ti) 



rWt2- r RM t,l3) 








Fig. 3. RM maps for all simulated patches ordered from high to low latitudes (top to bottom) in gnonomic projection. Their des- 
ignation and Galactic coordinates are quoted on top of each individual panel (also listed in Table 1). Different simulations for the 
same patch and for patches at the same latitude, but different longitudes, are shown in one row (see text for details). Each field has 
about 3000 x 3000 pixel and a size of about 1°.5 x l°5. The minimum and maximum of each map is shown together with the wedge 
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Table 3. Fitting results for structure functions and autocorrela- 
tion functions of RMs. 



name 



HN 

MN 

PNAC1 

PSAC 

MS 

HS 

PNC1 

PNC2 

PNAC2 



structure functions 



autocorrelation function 



AA 



Am 



34.8 
62.7 
2951.1 
2830.4 
64.9 
42.0 
13695.7 
31230.9 
2836.0 



0.4 
0.2 
29.9 
11.6 
0.5 
0.4 
293.8 
402.2 
57.4 



0.916 
0.862 
0.719 
0.751 
0.871 
0.909 
0.658 
0.676 
0.660 



0.006 
0.002 
0.007 
0.002 
0.005 
0.005 
0.010 
0.008 
0.013 



B AS 



As 



8.6 
15.7 
596.0 
596.8 
24.5 
11.5 
2255.7 
5792.6 
738.3 



0.1 
0.2 
3.9 
5.9 
0.4 
0.1 
16.4 
40.0 
5.5 



0.957 
0.853 
0.773 
0.831 
0.915 
0.913 
0.614 
0.661 
0.743 



0.008 
0.010 
0.013 
0.017 
0.010 
0.006 
0.018 
0.013 
0.014 



4.2.2. The structure function 

For turbulence studies, it is customary to investigate the power 
spectrum of the magnetic field, which can be retrieved from the 
power spectru m, autocorrelation fun c tion, or structure function 
of RMs (e.g. [ SimonettfetaD Il984t iMinter & Spanglerl 1 19961; 
lEnBlin & V ogt 2003). Observationally always uneven sampled 
RMs of EGSs are obtained. In this case the structure function 
is suited best (ISun & HanlliOOl) . Since RM maps are available 
as simulation results, the autocorrelation function is also investi- 
gated here as a crosscheck for the structure function results. The 
structure function (Drm(£#)) is defined as 



D RM (6e) 



[RM(0) - RM(0 + 6ff)Y 



(4) 



pairs 



where 60 is the angular separation of two pixels and A^ pa i rs is the 
number of RM pairs with the same angular separation. 

It would take a huge amount of computing time to obtain the 
structure functions from all pixels of our simulated maps. We 
simplified this by randomly selecting N s pixels to conduct the 
calculation for structure functions. These pixels are randomly 
distributed across the maps. To accurately assess the structure 
functions for small angular separations, the number N s should 
be as large as possible. As a compromise we selected N s = 2 x 
10 s pixels from the maps. The structure functions are shown in 
Figs . H El and fj (lower panels) . 

The theoretical deduction how structure functions are related 
to the power spectrum o f magnetic fields has been given by 
IMinter & Spangled dl996t) . where the structure function of RMs 
is found to follow a power law as Drm^S) = A6ff n with A being 
a constant and the spectral index m being related to the power 
low index a of P(k) as m = —a — 2. This relation holds for the 
inertial range, i.e., for the scale l\ < I < l a , where /; and l Q are 
inner and outer scales, respectively. The angular inertial scale 
can thus be defined as 50\ < 60 < 60 o , where 60{ = /;/r max and 
60 o = lo/r m ax- In principle the inertial range can be estimated 
by checking where the linear fit of the structure functions for a 
logarithmic scale fails. For the structure functions from our sim- 
ulated results, we are able to infer the outer scale but the inner 
scale is smaller than the pixel size. Below we name the outer an- 
gular scale 60 o as the transition angle. The fit results for structure 
fun ctions are listed in Tab le [3] 

EnBlin & Vogt (20031) have noted that the very inner part 
of the autocorrelation function of RMs (£rm) can be written as 
£rm(0) - £rm(80) = B60 s , where B is a constant, and s = -a - 2 
is the spectral index of the magnetic energy spectrum. This im- 
plies that the power law index of the structure function m and 
the power law index of the autocorrelation function s should 
be equal. We obtain the autocorrelation functions by a Fourier 



transform of the RM power spectra and show them in Fig. [5] The 
RM fit results for each patch are listed in Table[3] As can be seen, 
the spectral indices of the autocorrelation functions do not devi- 
ate very much from that of the structure functions. We account 
the differences mainly to computational effects when obtaining 
the autocorrelation functions. Unwanted components might be 
introduced during the pixelization of the maps when Fourier 
transformed. Therefore we are confident that the results based 
on the two methods are basically consistent with each other. 
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Fig. 4. PDFs and structure functions of the patches HN and MN 
at high and medium latitudes. See Tables Q] and [2] for details of 
the patch parameters. 



The spectral index of 5/3 is expected from the power law fit 
of the RM structure functions for a Kolmogorov-like turbulent 
magnetic field. This, however, is not seen in Table [3] Possible 
explanations are discussed in Sect. 4.7.1. 



4.3. Gaussian versus Kolmogorov-like random magnetic 
fields 



In the simulations presented by ISun et al .1 (I2008D . the random 
magnetic field components were realized in situ for each volume 
unit by Gaussian random numbers. The random fields created 
that way are called 'Gaussian fields' below. In that scheme, the 
simulated RMs are totally uncorrected between the map pixels, 
which implies that no small-scale structures are expected. When 
random magnetic fields with a Kolmogorov-like spectrum are 
generated, as described in Sect. 13.31 a large number of turbulent 
boxes are placed along the line-of-sight. For arcsec angular reso- 
lutions several neighbouring line-of-sights pass the same boxes, 
which results in extended coherent features in the simulated RM 
maps. 
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Fig. 5. Autocorrelation functions for RMs. 



For the northern patches HN and MN, we also performed 
simulations with Gaussian random magnetic field component for 
comparison. Both RM maps are shown in the first two rows in 
Fig. [3] and prove the above expectations. As seen from Fig. |4] 
the PDFs of the RM simulations with Gaussian random mag- 
netic fields nicely conform to Gaussian distributions, while the 
PDFs from the simulations with Kolmogorov-like random mag- 
netic fields do not. Quantitatively this can be verified by the 
skewness and kurtosis in Table [2] where the values of these two 
parameters are much larger for Kolmogorov-like magnetic field 
spectra than Gaussian ones. The structure functions shown in 
the bottom panel of Fig. |4]as well as the autocorrelation func- 
tions in Fig. [5] reinforce the argument. The structure functions 
for Gaussian random magnetic fields are flat, while the structure 
functions for a Kolmogorov-like magnetic field spectrum have 
slopes, which will be discussed below. 

4.4. RM patches from different realizations 

The RM structures in the present simulations are no predictions 
for sky emission structures towards the quoted directions, but are 
just one possible realization of the random magnetic fields fol- 
lowing a Kolmogorov-like power law spectrum in context with 
other model parameters as described before. 

We generate the random magnetic field components in the 
following way. Their amplitudes are determined according to 
the slope of the power-law spectrum with a spectral index of 
-11/3, and the phase angles are realized by random numbers. 
The rotation angle of the numerous 10 pc boxes placed along the 
line-of-sight are also determined from random numbers. These 
random numbers are actually pseudo random and are initiated 
by an integer called "seed". When the random "seed" number 
is changed, a different realization of the same random magnetic 
field distribution is obtained. 

To check structure variations and structure functions for dif- 
ferent realizations, we ran two more simulations for patch MS 



with different "seed" numbers. The RM maps are displayed in 
the fifth row of Fig. [3] and show different distributions of RM 
features as expected. Also the PDFs for the RM maps (Fig. |6]) 
vary considerably, as well as the parameters listed in Table [2] 
The structure functions of RMs, however, are almost identical 
as shown in the bottom panel of Fig.[6j except for regions with 
large angular separations exceeding about 10', where the struc- 
ture function becomes less accurate due to the decreasing num- 
ber of selected pixels. 

The mean RM for the three realizations of patch MS is be- 
tween -20 rad irT 2 and -28 rad irT 2 (Table |2| and the majority 
of RMs has negative values. However, due to the variance of up 
to 19 rad irT 2 (Table 13 all three simulations show a number of 
small extended features in the maps with positive sign. Naively 
this could be interpreted as reversals of the large-scale magnetic 
field in case the available number of observed RMs from EGSs 
is small. However, in contrast with large-s cale reve r sals i n the 
Galactic plane as discussed for instance bv lHan et al.l d2006l) . the 
reversals seen here are entirely caused by the distribution of ran- 
dom magnetic fields. For this medium-latitude area, we roughly 
estimate the probability for picking up a "reversal" to be about 
7%. This should be reflected in the distribution of observed RMs 
of EGSs and will prove the relevance of the simulations and the 
parameters used. 




69 (arcmin) 

Fig. 6. PDFs and structure functions of the patch MS for the 
same amplitudes but different realizations of the phases of the 
random magnetic field components created in a 10 pc box (see 
text for details). 



4.5. Galactic plane RM patches at different longitudes 

In this section we investigate RM variations in the plane versus 
Galactic longitudes. The four simulated patches have latitudes 
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of about +1°2. Two of them (PNC1 and PNC2) are in the inner 
Galaxy direction and the other two (PNAC1 and PNAC2) are in 
the anti-centre direction (see Table[T|i. 

The visible structures in the maps (the third row in Fig. [3) 
are less extended compared to those at higher latitudes. There is 
more fine structure for the inner Galaxy regions than for those in 
the anti-centre direction. The PDFs of all four simulated patches 
closely resemble Gaussian shapes (Fig. |7). The inner Galaxy 
patches PNC1 and PNC2 have smaller skewness values than the 
anti-centre patches (Tabled, indicating that their PDFs are more 
close to Gaussian shapes. Note that the skewness and kurtosis 
for all patches is very small, which indicates that it is impos- 
sible to infer whether the turbulence spectrum is Gaussian or 
Kolmogorov-like based on their PDF. The structure functions 
as discussed before, however, clearly reveal the properties of 
the underlying turbulence spectrum. The av erage RM value s are 
consistent with the RM results obtained bv lSun et all d2008l) . 

As displayed in Fig. [7] all Galactic plane structure functions 
are well fitted by power laws. Compared to the medium- and 
high-latitude fields a flattening of the structure functions occurs 
at smaller angular distances. The fit results for the power laws 
are listed in Table [3] From the parameters listed in Table [T] we 
obtain a transition angle of about 1' for PNC1 and PNC2 and 
about 3' for PNAC1 and PNAC2. Fig.[7]shows instead of a sharp 
change in the structure function a smooth transition from a linear 
slope to a flat curve in the logarithmic scale. The turn over of 
the structure functions ranges from about 1' to 2' for patches 
PNC1 and PNC2, from 3' to 5' for patch PNAC1, and from 2' 
to 5' for patch PNAC2, which are roughly consistent with the 
aforementioned expectations. The spectral indices are co nsistent 
with theoretical predictions by Minter & Spangler ( 1996) for 2D 
turbulence, which results in a spectral index of about 0.7 for the 
structure functions. However, our input turbulent magnetic fields 
are for 3D. The amplitudes of the structure functions increase 
towards the inner Galaxy regions, which is expected as the line- 
of-sights pass more turbulent cells or magnetic field boxes than 
towards the anti-centre or high-latitude direction. The amplitude 
of the structure function for PNC1 is larger than that for PNC2, 
which is mainly due to a larger electron density in this direction 
according to the NE2001 model. 



4.6. RM patches at different latitudes 

We finally compare the different results towards different lat- 
itudes. The patches HN/HS, MN/MS,and PNAC1/PSAC have 
longitudes between 130° and 140°. The latitudes are +70°, ±44° 
and +1°. For patch MS the "seedl" case is used. The six sim- 
ulated RM maps are all shown in the first column of Fig. [3] A 
clear trend for more extended structures for high-latitude regions 
is clearly seen. The PDFs for the high-latitude regions (HN/HS 
and MN/MS) all deviate significantly from Gaussian (Fig.[8]and 
Table©. 

The structure functions derived for the four high-latitude 
patches are shown in Fig. [8] and for the five plane patches in 
Fig. |7] All of them have similar properties. The transition an- 
gle is calculated to be about 10' for HN/HS and MN/MS re- 
gions, which is consistent with the results displayed in Fig. [8] 
The amplitudes of the structure functions are smaller for the 
high-latitude patches compared to those in the plane. This can 
be understood, because the length along the line-of-sightis small 
compared to that in the Galactic plane. The number of turbulent 
cells passed by the line-of-sights is smaller and hence is the am- 
plitude of the structure functions. We also note that the slopes 
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Fig. 7. PDFs and structure functions of the five simulated 
Galactic plane patches (see TableQ]for field parameters). 



of the structure functions increase for the high-latitude patches. 
These spectral changes will be discussed below. 

4.7. Discussion of RM simulations 

4.7.1. Theoretical structure functions revisited 

From the structure functions of simulated RM maps, we found 
evidence for latitude-dependent variations of their spectral slope, 
which gets steeper towards high latitudes. The regular magnetic 
field and the electron density have spatial irregularities, which 
are of larger scale than the sizes of the patches. Therefore the 
large-scale magnetic field cannot cause variations of the spectral 
slopes, but influences the amplitude of the structure functions. 
The coupling between th e electron density and random mag- 
netic fields introduced by ISun et al.l (|2008) augments the scat- 
tering of RMs and hence the amplitude of the structure function 
rather than changes its spectral index. For the present simula- 
tions, the random magnetic field is comparable to the regular 
magnetic field strength and dominates the regular field towards 
high-latitude regions. We conclude that all RM structures across 
the simulated patches are caused by the random magnetic field 
properties. 

The variations versus latitude must be ascribed to geomet- 
ric effects, which reflects in the integral length (rmax) as listed 
in Table [1] The size of the magnetic field box (L) is fixed in 
the simulations presented, therefore, the integral length deter- 
mines the number of turbulent cells or boxes passed by the line- 
of-sight. For a sufficient large number of cells, the fluctuations 
along the line-of-sight are averaged out. This results in a spectral 
slope of m — -a -3=2/3 for the structure function. Note that 
this spectral index was observed and interpreted as a 2D turbu- 
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Fig. 8. PDFs and structure functions of four patches (see first 
column of Fig. [3} at roughly the same Galactic longitude located 
at medium and high latitudes (see Table[TJ. 

Table 4. Fit parameters for structure functions from different 
simulation settings. 





L 


= 20 pc, 


''max = 10 kpc 




A foe) 




C 


AC 


m 


Am 


0.80 




3.304 


0.028 


0.727 


0.044 


0.40 




3.430 


0.021 


0.644 


0.033 


0.20 




3.387 


0.030 


0.760 


0.047 


0.10 




3.414 


0.009 


0.612 


0.014 


0.05 




3.447 


0.018 


0.622 


0.028 




A 


= 0.5 pc, 


'max = 10 kpc 




L(pc) 




C 


AC 


m 


Am 


200 




2.830 


0.006 


1.262 


0.006 


100 




3.045 


0.011 


1.081 


0.015 


50 




3.201 


0.008 


0.982 


0.010 


25 




3.401 


0.007 


0.709 


0.012 




A 


= 0.25 pc, L = 100 pc 




''max (kpc) 




C 


AC 


m 


Am 


0.01 




-6.270 


0.0181 


1.660 


0.014 


0.05 




-4.243 


0.0106 


1.413 


0.008 


0.10 




-2.689 


0.0128 


1.668 


0.009 


0.20 




-1.742 


0.0172 


1.355 


0.013 


0.40 




-1.206 


0.0202 


1.689 


0.013 


0.80 




0.368 


0.0133 


1.282 


0.012 


1.60 




1.630 


0.0088 


1.330 


0.009 


3.20 




2.711 


0.0060 


1.154 


0.006 


6.40 




2.970 


0.0083 


1.151 


0.010 


10.00 




3.083 


0.0078 


1.081 


0.010 



lence in the interstellar medium by iMinter & Spanglerl ([1996). 
2D turbulence is confined to thin sheets or filaments. However, 



the present simulations are based on a 3D distribution of the tur- 
bulent magnetic field components, which have not the shape of 
sheets or filaments. When the number of cells is small enough 
so that the integral length is comparable to the outer scale, the 
structure functi on tends to have the index of m = -a - 2 — 5/3 
as derived by (Mi nter & Spanglerl 1 19961) for 3D Kolmogorov- 
like turbulence. Therefore we conclude that in the inertial range 
(60\ < 59 < 60 o ) the derived structure functions shall behave like 
Z)rm(<56>) oc 6ff" with 2/3 < m < 5/3 depending on the number 
of cells. When the angular separation is larger than the transition 
angle, the structure function spectrum flattens and m gradually 
turns to zero. 

To prove these arguments we made several simulations for 
a patch of about 6° X 6° in the plane with A^side = 2048 corre- 
sponding to a resolution of about 1.'7. Although the resolution is 
much lower than the high-resolution simulations for the patches 
in Table [TJ the dependence of the structure functions should be 
the same. These simulations can be realized in much smaller 
computing time than the arcsec resolution simulations. We fit 
the structure functions by using log Drm(<S#) = C + mlogSO 
and list the results including errors of C and m in Table [4] We 
first fix the outer scale (L) and the integral length and vary the 
size of the small cubes (A). The results in Table [4] show that 
the size of the small cubes has little influence on the structure 
function parameter. We then change the outer scale, but keep the 
size of the small cubes and also the integral length constant. As 
expected the spectral index of the structure functions tends to 
be about 0.7 for very small outer scales, which means the num- 
ber of turbulent cells is large. Finally, we increase the integral 
length from 10 pc to 10 kpc by using the same values of cube 
size and outer scale. The spectral index of the structure func- 
tions varies between 1.3 and 1.7 until the integral length exceeds 
about 1600 pc, becoming smaller for larger lengths. These test 
simulations support our conclusion that the spectral index of the 
structure function primarily depends on the number of turbulent 
cells, which is mainly determined by the length of the line-of- 
sight. Thus, when assessing the turbulent magnetic field prop- 
erties from the structure functions of RMs, the line-of-sight de- 
pendence must be taken into account. In the Galactic plane the 
integral length is always large, which assures sufficient turbulent 
cells being passed through by the line-of-sight. Therefore the re- 
lation a = -m - 3 holds in general. Towards higher Galactic 
latitudes, however, the length of the line-of-sight decreases and 
hence also the number of turbulent cells decreases. The spectral 
index of the random fields approaches a = —m - 2. 

4.7.2. Diagnosis of the Galactic foreground 

It has been proposed to conduct a very dense RM survey of 
EGS with the SKA to study magnetic field structures in the 
Milky Way, in galaxie s, in cluster of galaxies and beyond 
dBeck & Gaenslerl2004l) . The key issue to achieve relevant infor- 
mation for extragalactic objects is to properly extract and sepa- 
rate the foreground RM contributions from the observations. We 
have tested the minimum number of pixel needed to calculate 
the correct structure function for the simulated maps and found 
that about 10000 pixel are required. For our map size this corre- 
sponds to a mean angular separation of about V in case the EGS 
are not larger than the pixel size of l'/6 and their intrinsic RMs 
average out. 

In the plane the PDFs are close to Gaussian. Therefore the 
average of RMs is a good estimate of the foreground contribu- 
tion caused by the large-scale ma gnetic field and t he thermal 
electron distribution as modelled by lSun et al.l ((2008). However, 
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large variances (Table [2]i and amplitudes of the structure func- 
tions (TableO imply that after having subtracted the average RM 
in a certain direction the RM scattering of the foreground emis- 
sion is still preserved. Thus patches in the plane are certainly not 
well suited for the study of extragalactic magnetic fields. 

At high latitudes the PDFs deviate significantly from a 
Gaussian distribution and also more extended structures show 
up (Fig. |3). To take the foreground RMs appropriately into ac- 
count a dense sample of EGS is required. In case an extended 
extragalactic object is seen towards an excessive RM patch 
of the same size, additional information is required to solve 
the RM foreground problem. On the other hand, the variation 
of the RMs are much smaller than in the plane. Thus high- 
latitude regions seem better suited to study extragalactic mag - 
netic fields. However, according to Kronberg et al. (2008) cos- 
mological magnetic fields of a few 10~ 7 /uG at a redshift of 2 have 
an expected RM excess of up to 20 rad rrT 2 . This is of the same 
order of magnitude as the RM variance of the Galaxy at high 
latitudes obtained from the present simulations, which means a 
detection of cosmological magnetic fields is quite challenging. 
Note that reducing the coupling fa ctor between ele ctron den- 
sity and random fields introduced bv lSun et al.l (|2008) to explain 
the observed depolarization of Galactic synchrotron emission at 
low angular resolution will in turn also reduce the foreground 
RM variance. Therefore future high-resolution observations of 
the small-scale structures of the magnetized ISM are needed to 
investigate the coupling factor in detail. Such observations are 
planned to be carried out with LOFAR, which will operate at 
m-wavelengths and will thus be able to detect very small RM- 
variations from small-scale Galactic emission features. 

5. Simulated total intensity and polarization maps 

We briefly describe the properties of the simulated total inten- 
sity and linear polarization maps emerging from the same sim- 
ulation runs where the already discussed RM distributions were 
calculated. These simulated Stokes /, U, Q maps are all for the 
frequency of 1 .4 GHz, whereas the RM distribution is frequency- 
independent. 

The simulated total intensity /, U, Q, polarized intensity PI 
and polarization angle PA maps are shown in Figs. |B.l||B.2||B.3l 
IB. 41 and IB.5l for the same patches and simulation variations as 
displayed in Fig. [3] for RMs. The statistical parameters for the 
/, PI, and PA maps are calculated and listed in Table [5] The 
structure functions for / and PI are shown in Fig. [9] 

Similar as for RM maps, the patches in the plane are dom- 
inated by structures on very small scales, while the patches at 
high latitudes exhibit more extended features. These structures 
are all related to the properties of the turbulent magnetic fields. 
The variation of the sizes of the features is primarily due to the 
much larger integral length along the line-of-sights in the plane 
than at high latitudes. 

The large-scale cosmic -ray and magnetic field also have im- 
prints on the simulations. For example, the clear gradient of the 
total intensity particularly in the plane maps towards the inner 
Galaxy traces the large-scale magnetic field and the cosmic -ray 
distribution. The asymmetric halo field adds up to the disk field 
and enhances the field strength below the Galactic plane, which 
produces larger intensities there. This can be seen from the av- 
erages of / and PI intensities listed in Table|5]for the PSAC and 
the PNAC1 patch. 

The average polarization percentage is about 10% for the HN 
and MN patches and about 16% for the HS and MS patches. 
Towards these directions both RM and its variance are small 



so tha t Faraday depolarization is small as well (Sokoloff et al. 
119981) . The low degree of polarization is thus due to the ran- 
dom magnetic fields, which dominate the field components. The 
percentage becomes about 7% for patch PNAC1, 5% for patch 
PNAC2, and 2% for the inner Galaxy patches, which is caused 
by strong Faraday depolarization, as both the RM and the RM 
scattering are large. There is no indication for a correlation 
among the variances of RM, polarized intensity PI and polar- 
ization angle PA. 

The structure functions for total intensity / (Fig. |9]l have a 
linear slope in the logarithmic scale, but these are generally more 
shallow than those for the corresponding RMs. The total inten- 
sity / and polarized intensity PI maps at high latitudes have quite 
similar structure functions, while they are rather distinct from 
each other in the plane. The structure functions of PI are nearly 
flat. This indicates that due to strong depolarization in the plane 
the PI structure functions do not contain the fluctuation infor- 
mation of the turbulent magnetic field any more. 

6. Conclusions 

Galactic 3D-models of the distribution of thermal electrons, 
cosmic -ray electrons, and magnetic fields were used to calculate 
total intensity, polarization and RM maps for selected patches of 
the Galaxy via the hammurabi code. The angular resolution of 
the maps is about 1'.'6, close to that of the planned SKA. The 
simulated patches are distributed at different longitudes and lat- 
itudes, which enables us to study spatial variations of the re- 
sults. The simulated RM maps are of particular interest, as they 
are directly related to the planned deep RM EGS survey with 
the SKA. Complete SKA simulations need in addition simulated 
sky patches wi th distributed polariz ed sources. Such simulations 
were made bv lWilman et alj d2008l) . 

In the present simulations the random magnetic field com- 
ponent is assumed to follow a Kolmogorov-like power spectrum. 
Compared to purely random Gaussian magnetic fields, structures 
on different scales show up in the RM as well as in the total in- 
tensity / and the polarized intensity PI maps. The simulations 
predict more extended structures for high-latitude regions. 

We study the PDFs and structure functions for the RM maps. 
The PDFs of the patches in the plane are close to Gaussian, 
whereas at high latitudes clear deviations from Gaussian shapes 
are noted. The structure functions have a power law slope for 
angular separations smaller than the transition angle and flat- 
ten for larger angular separations. The transition angle is deter- 
mined by the outer scale of the turbulence and by the length of 
the line-of-sight. The amplitudes of the structure functions be- 
come smaller and the slopes of the structure functions steepen 
towards higher-latitude regions. All this is entirely caused by ge- 
ometric effects. In the plane, the integral length is much larger 
than at high latitudes, therefore more turbulent cells are passed 
by the line-of-sights. This significantly smears out fluctuations 
along the line-of-sights, which makes the 3D turbulence obser- 
vationally a 2D turbulence. Thus the structure functions are flat- 
ter in the plane than at high latitudes. We notice that the PDFs 
of RMs have a wide distribution in all directions, which poten- 
tially leads to depolarization of Galactic radio emission, when 
observed with larger beams. It also strongly influences the in- 
terpretation of discrete features sampled by RMs of EGS. It is 
obvious that only a very dense grid of extragalactic RMs will 
help to separate Galactic foreground influence from extragalac- 
tic contributions. This is planned for the SKA. 

The simulations at 1 .4 GHz presented in this paper can be 
easily run at any other frequency. While the RM maps are fre- 
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Table 5. Statistical parameters for the simulated /, PI, and PA patches. The patch designations are listed in Table [TJ 
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variance 
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(mK) 
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0.31 


-0.049 


48.59 


13.63 


-0.04 


-0.347 


32.85 


8.81 


-0.36 


0.191 


HN (Gaussian) 


446.73 


13.98 


0.06 


0.003 


57.91 
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0.03 


0.003 
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6.15 
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0.122 
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489.76 
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0.004 
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2.527 


PSAC 


1630.36 


39.12 


-0.28 


0.233 


111.65 


50.87 


0.19 


-0.554 


-18.93 


23.68 


0.54 


2.327 


MS (seedl) 


519.56 


18.59 


0.46 


-0.079 


85.53 


25.60 


-0.78 


-0.101 


-35.46 


18.62 


0.82 


4.854 


MS (seed2) 


540.35 


16.87 


0.18 


0.873 


122.55 


25.25 


-0.65 


0.170 


-29.04 


14.70 


-0.02 


-0.686 


MS (seed3) 


556.77 


19.11 


-0.09 


-0.563 


87.20 


19.04 


-0.12 


-0.296 


-37.31 


11.56 


0.39 


-0.202 


HS 


512.79 


24.02 


0.25 


-0.492 


78.68 


22.62 


-0.15 


-0.290 


-24.73 


21.44 


0.04 


0.336 


PNC1 


2996.47 


121.62 


0.10 


-0.710 


61.85 


32.01 


0.59 


0.140 


-14.66 


50.31 


0.54 


-0.850 


PNC2 


2780.71 


93.90 


0.13 


-0.040 


54.90 


27.98 


0.56 


0.075 


15.26 


37.16 


-0.66 


0.348 


PNAC2 


1449.76 


40.76 


0.31 


0.419 


71.01 


29.17 


0.19 


-0.235 


45.83 


31.16 


-2.51 


7.768 



* o 

t o 



: . 4r x 



PNAC2 
PNC2 
PNC1 
PNAC1 



100 



1 

68 (arcmin) 




(arcmin) 



*********** * * * ;f *.$ 4** 



° ° 9 S x x 



q J Q 8 x> $ ' 



******* 



PNAC2 
PNC2 
PNC1 
PNAC1 



* * * * 6 
666« 



100 



10 



ioooooooo 



ooo 



x x x ^ 
ooooo 



o0 oooooo° 
□ DD n 

n . A A AqAq 



□ 



□ 



□ 



n ° " 
□ A 

□ a i i* 

□ A _L * 

□ A a_ W 

□ A + * 



HS 

MS 

PSAC 

PNAC1 

MN 

HN 



1 

68 (arcmin) 



0.1 



10 



68 (arcmin) 



Fig. 9. Structure functions for the total intensity / (top panels) and the polarized intensity PI maps (bottom panels). 



quency independent, total intensity and polarization simulations 
are of particular interest for low-frequency simulations, where 
LOFAR is expected to do sensitive high-resolution observations 
in the near future. At low frequencies, thermal absorption adds to 
depolarization effects and must be properly taken into account. 
LOFAR will make a major attempt to trace the very weak signal 
from the Epoch of reionization, where the selection of a high- 
latitude Galactic area with very low foreground signal is impor- 
tant for the success of these ambitious low-frequency observa- 
tions as well as the subtr action of a l l kind s of extragalactic fore- 
grounds as discussed bv lJelic et alj d2008l) . 
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Appendix A: Details for the realization of random 
magnetic field components 

The general concept is to generate random magnetic fields in 
a box by a Fast Fourier Transformation (FFT), randomly rotate 
the box around its centre to obtain a different view on the box, 
and finally fill the volume of the Galaxy with these randomly 
rotated boxes. The size of a box is L 3 and the size of the Galaxy 
is D x xD Y x D Z - 

We first define several Cartesian coordinate systems: 

- (X, Y, Z): with XT-plane parallel to the Galactic plane and 
the Galactic centre located at (Dx/2, Dy/2, Dz/2). 

- (x, y, z): based on the box with the origin at its vertex. When 
the Galaxy is filled with boxes the origin of the first box is the 
same as for the system (X, Y, Z) and the xy-plane is parallel 
to the .XT-plane. 

- (x', y' , z') or (x", y", z"). the same as (x, y, z), but with the 
origin at the centre of the box. 

Based on these coordinate systems the whole process is il- 
lustrated in Fig. [AT] 

We first consider the random magnetic field distribution in a 
box. The random magnetic field is assumed to be isotropic and 
homogeneous, and its three components do not correlate with 
each other. The three components thus follow the same power 
law spectrum P(k) - C^k" . Here is the amplitude, k - 2n/l\s 
the wave vector with / being the spatial scale and a is the spectral 
index. For the simulations presented in this paper a Kolmogorov- 
like turbulence is assumed, which means a = -11/3. 



We first generate the ID component of the random magnetic 
field. The other two components were obtained by simply re- 
peating the process for the first one. This ID component denoted 
as b(x, y, z) is the inverse Fourier transformation of the power 
spectrum as the following, 



b(r) 



--f^f(PW) 



exp(/£(r) + /k • r) dr, 



(A.l) 



where £ is the phase angle taken as a uniform random number 
between and 2n, i = V— T and r = (x, y, z) is the position vec- 
tor. In practice Equation ( I A. It needs to be discretised. Therefore 
we divide the box into A' 3 cubes with a size of A 3 with A = L/N 
(Fig. IA.lt . Equation (IA. lb is then expressed as, 



N-l N-l N-l 

h'm'n' = 2^ Zj XI ' Pl " m X eX P( /A ^' / ' + k >" m ' + k n n '))> 
1=0 m=0 /?=() 



(A.2) 



where *P, 



him 



exp(/£ m „), P 

r 



l/iUl 



4 



h = 



tf+kl+ k 2 „. Here 



2nl/L 

2n(N - 1)1 L 



for I < N/2 
for I > N/2, 



= C l k im„> and k i™ = 



(A.3) 



with the same for k m and k„. Generally ki, n „ is limited to the in- 
ertial range 2n/l < ki mn < 2n/h, where k — A and l a = L are 
the inner scale and outer scale of the turbulence, respectively. To 
assure that by m i n ' is a real number r Pi„ ln should be Hermitian, i.e. 
conforming to the following, 



- 'Plmn - PN-lN-mN-n, 

- 'Plmn is a real number, if the following is satisfied 



I — m — n = 

/ = 0, N/2, m = 0, N/2, n = 0, N/2 



(A.4) 
mod(M 2) = 1 



mod(7Y, 2) = 0. 



(A. 



The FFT package provided by the FFTW library is applied 
to obtain the random magnetic field distribution. 

To calculate the random magnetic field in any position (X, Y, 
Z) in the Galaxy b(X, Y, Z), the following procedure is executed 
(Fig. ED, 



1. transform (X, Y, Z) to (x', y', z') as, 

. X L Y L , Z L ^ 

x — i v L , y = i v L , z ( 7 L (A. 6) 

L 2 L 2 L 2 ' 

where i x = mtiX/L), i y = int(F/L) and i z = int(Z/L) are the 
integer part. 

2. rotate (x', /, z') by Euler angles {if/, 8, <p) around the centre 
of the box to get (x", y" , z") as, 



(A.7) 



x" ' 




x' 


y" 


= BCD 


y 


_z" _ 




. z' . 



and 



B = 



cos(/( sini/r 
- sin i// cos x// 
1 



(A.8) 



www.rFtw.org 



14 



X. H. Sun & W. Reich: High resolution simulations 





y=jt,e=o,it>=o 































































i|/=37t/2,e=0,(|>=0 





























































Appendix B: Simulated total intensity and 
polarization maps 

We display the total intensity / maps in Fig. IB. II Stokes U maps 
in Fig. IB. 21 Stokes Q maps in Fig. IB. 31 polarized intensity PI 
maps in Fig. IB.4l and polarization angle PA maps in Fig. 
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Fig.A.l. 2D-scheme for the realization of random magnetic 
fields. The sizes of the box and the cube are L and A, respec- 
tively. These boxes cover a region of Dx x Dy- The Euler angles, 
by which the boxes are rotated from the lower left box, are indi- 
cated. The shaded box marks the same cube after rotation. 



1 
C = cosO sin6» 
- sin 9 cos 6 

cos<p sin^i 
D = - sin <p cos <p 
1 
Here the Euler angles are obtained as, 



(A.9) 



(A. 10) 



and 
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n 

2 

f o 
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n 
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£ 

2 
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2£ 

2 



0< 7?i<0.25 
0.25<-Ri<0.5 
0.5<7?i<0.75 ' 
0.75<"Ri<l 

0< "R 2 <0.33 
0.33< n 2 <0.67 , 
0.67<7? 2 <1 

0< ??3<0.25 
0.25< ??3<0.5 

0.5<7? 3 <0.75 
0.75<??3<1 



(A. 11) 



(A. 12) 



(A.13) 



Here we use three random numbers %2, and to deter- 
mine the Euler angles. For simplicity the Euler angles are set 
to be integer multiple of n/2 to avoid complex projections. 
3. transform (x", y", z") to (x, y, z) then to (/', m', n'), 



L 

2' 



y = y + 



r 



■ z + 



2' 



and 



V = int(-), 



int(l), 



int(-). 

A 



(A. 14) 



(A.15) 



Finally b(X, Y, Z) = by, 
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Fig. B.l. Total intensity / maps for all the simulated patches in gnomonic projection. The centre coordinates are given on top of 
each panel together with its label (see Table[T]for details). Intensities are in mK brightness temperature at 1.4 GHz. The minimum 
and maximum of each man is shown together with the wedee below each nanel. 




Fig. B.2. As Fig. IB. 11 but for the simulated Stokes U maps. 
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Fig. B.3. As Fig. IB. II but for the simulated Stokes Q maps. 






Fig. B.4. As Fig. IB. 11 but for PI maps calculated from the simulated U and Q maps. 
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Fig. B.5. As Fig. IB. II but for the polarization angles PA calculated from the simulated U and Q maps. 



